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Abstract 

Statistical and electrical properties of the conduction electrons of a silver or gold sphere with a radius from 1 to 2 nm 
are shown to differ drastically from the properties of electrons in a bulk metal sample. If the radius of a noble metal 
sphere decreases from 1 0 to 1 nm, its conductivity oscillates around the bulk metal value with increasing amplitude 
and drops at the 'magic' numbers of electrons. These numbers are equal to 1 86, 1 98, 254, 338, 440, 556, 676, 832, 91 2, 
1 ,284, 1 ,502, and 1 ,760, in agreement with various experimental data. We show that the conductivity and capacitance 
of a metal nanosphere can be changed by several orders of magnitude by adding or removing just a few electrons. 
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Background 

In the past, a measurement of optical absorption by silver 
nanoparticles embedded in glass showed that the particles 
had normal metallic properties when their diameters were 
decreased down to 2.2 nm [1], Contrary to this finding, 
metal particles with sizes below 2 nm cannot be conduct- 
ing according to more recent papers [2,3]. Very recently, it 
was understood that the metal-insulator transition (MIT) 
is gradual so that particles with certain magic' numbers 
of electrons become insulating while others remain con- 
ducting [4]. If electrons move inside a sphere, then the 
numbers 186, 198, 254, 338, 440, 556, 676, 832, 912, 1,284, 
1,502, and 1,760 are known to be magic'. It was experi- 
mentally found that the above numbers are indeed magic 
for clusters of many metals [5-16]. This allows one to 
consider the motion of electrons in a spherical jellium 
[8,12,17,18]. 

We recently studied statistical properties of 500 to 2,000 
free electrons confined in a spherical potential well with a 
radius from 1.2 to 2 nm. The averaged occupation num- 
bers of the electron energy levels and the variances of 
the occupation numbers were computed for both iso- 
lated metal nanoparticles and those in equilibrium with 
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an electron bath. The sum of the variances of all occu- 
pation numbers was found to depend on the number of 
electrons nonmonotonically dropping by a few orders of 
magnitude at 'magic numbers' of electrons. Here, we show 
how the statistical properties of the conduction electrons 
are related with the electrical properties of metal nanopar- 
ticles. Calculations of the DC conductivity and capaci- 
tance of single nanometer-sized noble metal spheres are 
reported. We predict a transistor-like behavior of a sin- 
gle nanoparticle when an additional charge of the particle 
drastically changes its conductivity and capacitance. 

Methods 

Statistical and transport models 

The electron statistics and capacitance of metal nanopar- 
ticles are investigated by the Gibbs ensemble method. 
The transport properties of electrons are studied using a 
quantum-mechanic generalization of a formula that can 
be obtained by solving the Boltzmann transport equation. 

The Boltzmann transport equation (BTE) is widely used 
in modern physical kinetics [19,20] even at nanoscales 
[21,22], In our previous work, we solved BTE and deter- 
mined the response of the conduction electrons on an 
electromagnetic wave allowing for dependence of the dis- 
tribution function on the wavenumber. Then, we obtained 
the spatially dispersive permittivity of metal and applied 
a generalized Mie theory [23] to calculate spectra of light 
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extinction by silver nanospheres. Because of excitation of 
the rotationless (longitudinal) plasmon-polariton waves, 
the frequency of the Frohlich resonance co = 3.5 eV [24] 
was found to increase with decreasing the radius a, so co = 
3.635 eV for particles with a = 1.5 nm and co = 3.73 eV 
at a — 1 nm (see Figure two of [25]; a similar depen- 
dence of co on a as well as the formulas of [23] is found 
in a recent paper [26]). The blueshift by 0.15 eV and the 
width of the plasmon resonance calculated for a particle 
beam created by Hilger, Tenfelde, and Kreibig [27] were in 
excellent agreement with the experimental data. We con- 
cluded that silver clusters with a < 1 nm do not contribute 
into the extinction spectra of the beams. However, it was 
not possible to establish whether the contribution of these 
ultrathin particles into the integral extinction spectrum 
vanished due to MIT. 

Energies of the conduction electrons 

It is a common practice to assume that the conduction 
electron energy e in metal is a function of the continuous 
quasi-momentum p = fik. This approach can be used to 
model the properties of bulk metals in which an electron 
state s can be described by a set of four quantum num- 
bers, s = s z , p, where s z is a projection of the electron 
spin. However, if an electron moves inside a microscopic 
sphere, the vector p can no longer describe an electron 
state. The set p x >Py> and p z has to be replaced by a set of 
the radial q, angular momentum /, and magnetic m quan- 
tum numbers. Then, integrals over p should be replaced 
by sums over the electron states s, according to the follow- 
ing rule: J2 S = ^ f dp, where V is the volume of the body 
and h is the Plank constant. In this study, we used discrete 
sets of the electron energies e s . 

Shape of metal nanoparticles 

Every set of the electron energies e s was obtained at a cer- 
tain number N of electrons confined in a potential well of 
an ideal spherical shape. We used the parameters of bulk 
silver and gold [28], namely, the depth of the potential well 
U = 9.8 eV and the Wigner-Seitz radius r s = 0.16 nm. 
Because of the spherical symmetry of the electronic wave 
functions, different physical processes have peculiar fea- 
tures at magic numbers of electrons. The magic numbers, 
which we determined for perfect noble metal spheres, 
are presented in the 'Background' Section of this paper. 
Surprisingly, almost the same values were determined in 
various experiments working with particles of different 
metals, even though no special efforts were made to pro- 
duce perfect spheres. In addition, different theoretical 
papers also reported similar magic numbers, according to 
Figure 1. This means that effects associated with the pecu- 
liarities of the spacing of e s in spherical nanoparticles are 
sensitive neither to surface distortions nor the values of 
the parameters U and r s . 



Results and discussion 

Variances of the occupation numbers 

In our previous work [29], we reported statistical prop- 
erties of the conduction electrons in isolated metal 
nanospheres. To study the systems with a fixed number 
of electrons, the method of the canonical ensemble was 
applied. The averaged occupation numbers (n s ), variances 
of the occupation numbers {n 2 } — {n s } 2 , and sums of the 
variances A = J2 S ((n 2 ) — (n s ) 2 ) were computed and dis- 
cussed. In [29], we also examined the properties of the 
conduction electrons in grand canonical ensembles where 
the chemical potential fio was fixed. Figure 2 represents 
the values of A calculated at fixed N (canonical ensem- 
bles) and /jLo (grand canonical ensembles). The sum of the 
variances depends on the number of electrons nonmono- 
tonically dropping by several orders of magnitude at magic 
numbers of electrons. The decrease in A can occur if (i) 
the distance between the Fermi level and the neighboring 
higher energy level, ey + i — €f, is large compared to the 
thermal energy and (ii) the Fermi level is fully occupied at 
absolute zero temperature. Addition of one atom to a par- 
ticle with N m conduction electrons results in a substantial 
increase in the Fermi energy, as is evident from Figure 2a. 
If a particle has a magic number of electrons, the chemical 
potential lies in the gap between the distant energy levels, 
so the number of the current carriers is greatly reduced. 
The influence of this effect on the electrical properties of 
the metal nanoparticles is studied below. 

Conductivity 

The response of the conduction electrons of metals to 
an infrared and far infrared radiation is well described 
by a Drude dielectric function [30]. In the corresponding 
limit of small emission wavenumbers, this function can be 
derived by using either a quantum theory by Lindhard [31] 
or the classical Boltzmann transport equation [32] (see 
derivations in [20]). The permittivity of metal is expressed 
through the conductivity a that, in the DC case, can be 
calculated from the following equation: 

o = --—n e J , (1) 

where e and m e are the charge and effective mass of the 
electron, respectively, n e = N/V, e = p 2 /(2rn e ),f(e) is an 
equilibrium electron energy distribution function, and r 
is the relaxation time. From Equation 1, the classical result 

2 

a b = ^- n e x is obtained at r = const and any/(6) finite 
at 6 = 0 and vanishing at e —> oo. The formula for <Jb can 
also be derived by substituting a zero-temperature Fermi- 
Dirac distribution function into Equation 1. 
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Figure 1 Experimental (centered boxes with error bars) and theoretical (crosses) 'magic 7 numbers of electrons in metal clusters. Solid grid 
lines indicate A/ m = 1 86, 1 98, 254, 338, 440, 556, 676, 760, 832, 91 2, 1 ,01 2, 1 ,1 00, 1 ,284, 1 ,502, and 1 ,760. Dashed grid lines indicate A/ m = 268, 542, 
1,074, and 1,206. 



A generalization of Equation 1 for discrete energy levels 
gives the following formula: 



2 

3 m P 



V ^ 



kT 



[1 - (If,)] (If,) T (€,) , 



(2) 



where (n s ) is the averaged occupation number of the state 
s. We tested Equation 2 by computing the normalized 
conductivity defined at constant r, 



6(A0 = ^ = 3^£^ [1 - <Ws>] {Hs) - 



Oh 



kT 
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The equality 6 = 1 should hold for large' particles since 
properties of a macroscopic body are independent of the 
boundary conditions for the electron wave function. The 
calculations were performed by using sets of e s for N free 
electrons confined in a spherical potential well with the 
radius a = r s A/ 4 / 3 , where r s = 0.16 nm. Figure 3a presents 
the results obtained at N in the range from 2,000 to 2.5 x 
10 5 , T = 300 K. There are pulsations of (3 vanishing as 
sphere radii increase above 9 nm that corresponds to N > 
2 x 10 5 . Therefore, Equation 2 works well, and particles 
with a > 10 nm can be regarded as macroscopic. The left- 
hand side of the curve in Figure 3a (at a from 2 to 4.5 nm, 
i.e., N from 2,000 to 20,000) shows the oscillations of (5(a) 
with the amplitude increasing with the decrease of a. 

The conductivity 6 at Af = 200 to 2,000 was calculated 
by using more realistic values of e s found for a spheri- 
cal potential well with the parameters of silver and gold. 
According to Figure 3b,c, the value of (3 is not a mono- 
tonic function of N and drops sharply when N is equal to 
one of the magic numbers Af m . The appearance of magic 



numbers is a general property of fermionic systems. In this 
paper, the magic numbers of the conduction electrons are 
identified by the dips in the conductivity <5(N). The values 
of N m and 6 (N m ) x 100 are listed in Table 1. The found 
values of N m are in excellent agreement with the experi- 
mental and theoretical magic numbers of clusters of many 
metals according to Figure 1. 

All the experimental numbers N m in Figure 1 were 
obtained by using the mass spectroscopy from dips in the 
mass spectra. For example, Katakuse and co-workers [6] 
found magic numbers of atoms equal to 197 for negative 
cluster ions of silver (Ag)~ and 199 for positive (Ag)+ 
cluster ions. Other magic numbers of atoms were 137 
for (Ag)-, (Au)-, and (Cu)" and 139 for (Ag)+ (Au)+ 
and (Cu)+. In these cases, the negatively and positively 
charged and neutral particles had the same magic num- 
bers of electrons, N m = 198 and 138. Thus, the anomalous 
properties of the metal nanoparticles in the experiments 
[5-15] are determined by electron motion [29] but not 
their atomic structure. Moreover, the model of single elec- 
trons trapped in a spherical potential well was shown to 
be adequate [6] though the shape of the clusters obtained 
by the bombardment of metal sheets with Xe ions was not 
controlled. 

A nonlinear dependence of 6 on N can occur even in a 
single sphere if N varies around Af m . To examine electric 
properties of a single charged nanoparticle, let us consider 
a sphere in thermal equilibrium with a reservoir of elec- 
trons, so the electrochemical potential \i = /xn + e<p is 
constant inside the sphere; here /xn is the chemical poten- 
tial of the neutral sphere and 0 is the electric potential. 
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Figure 2 Fermi energies and variances of the occupation numbers 
of electronic states of single Ag or Au spheres, (a) Fermi energy as 
a function of the number N of conduction electrons, (b) Sums of the 
variances A normalized to the bulk metal value Ab in canonical 
ensembles (points) and grand canonical ones (crosses). The grid lines 
are the same as in Figure 1 . 



For a fixed /x, we determined N = ^2 s (n s ) by using 
the Fermi-Dirac occupation numbers and computed the 
charge of the sphere Q = e (N — Nq) , where No is the 
number of electrons in the neutral particle. We calculated 
the quantities Q and 6 for a charged 336-atom Ag or Au 
nanoparticle. We found that the 336-atom particle holds 
two extra electrons when the value <p changes in a wide 
range of about 0.6 V. If the mean number of electrons in 
the particle is equal to 338, then 6 = 4.4 x 10 -5 . The 
normalized conductivity of the neutral sphere is found to 
be 6 = 0.78; 6(338)/6(336) = 5.7 x 1(T 5 . In the con- 
sidered example, the neutral sphere is conductive, but the 
charged one with two extra electrons turns out to be an 
insulator. 

Capacitance 

A parameter that describes the dependence of Q on 0 is 
the electric capacitance 



A straightforward calculation of the derivative of Q gives 
the capacitance of the charged particle with 338 electrons 
C = 6.1 x 10 -22 F that is much lower than C = 1.1 x 
10 -17 F of the neutral 336-atom sphere. The change in the 
capacitance C(338) / C(336) = 5.3 x 10 -5 is similar to the 
the correspondent change in the conductivity. 
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Figure 3 Normalized DC conductivity, (a) Normalized DC conductivity vs rigid-wall sphere radius a = r s A/ 1/3 at N from 2,000 to 2.5 x 10 5 . 
Normalized DC conductivity of a neutral silver or gold sphere at (b) N = 180 to 382 and (c) N = 382 to 2,000. The grid lines are the same as in Figure 1. 
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Table 1 Normalized conductivity (%) calculated for an Ag 
or Au particle with a magic number of atoms 









N m 












186 


198 


254 


338 


440 


676 


912 


6 (%) 


0.03 


2.6 


0.01 


0.005 


0.37 


5.6 


4.1 



By calculating the derivative of Q in Equation 4 at A/" 
defined through the Fermi-Dirac occupation numbers, 
we get 



where A is the sum of the variances of the occupa- 
tion numbers shown in Figure 2 by crosses. Equation 5 
expresses the relation between the reaction of the con- 
duction electrons to the electric field and the fluctuations 
of the occupation numbers of the electron states. Thus, 
the peculiarities of spacing and degeneracy of the elec- 
tronic energy levels have similar effects on the statistical 
and electrical properties of a nanometer-sized particle. 

During the calculations we neglected Coulomb effects. 
These effects are as follows. When an electron leaves 
a neutral metal sphere, it overcomes the attraction of 
the positive charge remaining on the sphere. Conse- 
quently, the work function increases by the value AU = 
0.54/^(nm) eV [33]. For example, AU ~ 0.5 eV for a 
338-atom noble-metal sphere. In addition, Coulomb 
forces should change the binding energy of the electron 
in a charged sphere compared to that in the neutral one. 
Though the change in U can be large, it should not be 
critical to the effects studied in this paper. Indeed, they 
depend on the value of 6y + i — €f, but this difference is a 
weak function of U For example, for a noble metal sphere 
with 338 conduction electrons, we get ey + i — €f = 0.69 eV 
at U = 9.8 eV, and e /+ i - e f = 0.74 eV if U -> oo. 

Conclusion 

In conclusion, the statistical properties, conductivity, and 
capacitance of a single nanometer-sized metal sphere 
depends very strongly on the number of conduction 
electrons N in the range from 200 to 2,000. In par- 
ticular, the DC conductivity drops by several orders of 
magnitude if N is equal to one of the magic numbers. 
For instance, addition of two electrons to a 336-atom 
noble metal sphere should reduce both the conductiv- 
ity and capacitance of the particle by four orders of 
magnitude. 
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